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Abstract 

The hypothesis of a tiny fraction of the cosmic inventory evolving cosmologically as a degenerate 
Fermi gas test fluid at some dominant cosmological background is investigated. Our analytical 
results allow for performing preliminary computations to the evolution of perturbations for rela- 
tivistic and non-relativistic test fluids. The density fluctuation, 5, the fluid velocity divergence, 8, 
and an explicit expression for the dynamics of the shear stress, a, are obtained for a degenerate 
Fermi gas in the background regime of radiation. Extensions to the dominance of matter and to 
the ACDM cosmological background are also investigated and lessons concerning the formation of 
large structures of degenerate Fermi g depicted. 
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I. INTRODUCTION 



The complete understanding of the cosmic inventory is one of the most challenging issues 
of theoretical physics. The efforts to relate theory with the observational data and to 
achieve an enlarged overview of the dark sector do necessarily require an interplay between 
general relativity, astrophysics and particle physics [l-4]. It may provide guidelines to the 
resolution of some fundamental questions that indeed influence the development of particle 
physics beyond its standard model (SM). 

Our proposal concerns the investigation of matter components evolving cosmologically 
as a degenerate Fermi gas (DFG). We apply the linear cosmological perturbation theory to 
obtain an analytical description of the growth rate of perturbations for a DFG test fluid. It 
is suggested that particles behaving as a DFG are related to a tiny fraction of dark matter 
(in fact, less than 1%) that could be made of massive neutrinos in adiabatic equilibrium 
with the radiation background Q|. In fact, just under very particular circumstances it 
provides attractive conditions for a model of dark matter components [3, Isj]. In general, 
neutrinos in the cosmic background are treated either in the limit of radiation regime or 
in the limit of cold matter. Treating massive neutrinos as a DFG allows for quantifying a 
smooth transition between ultra-relativistic (UR) and nonrelativistic (NR) thermal regimes. 
It becomes a convenient tool used for describing the role of such particles in the cosmic 
inventory, from radiation-dominated (RD) to matter-dominated (MD) eras. 

In the context of the standard cosmology, theoretical speculations about the existence of 
a degenerate Fermi gas (DFG) of neutrinos produced in the early universe were carried out 
in the past 9|. Depending on the value of the mass, the density of neutrinos is enhanced 
when they cluster into gravitational potential wells, although that overdensity is limited 
to small factors [10||. Otherwise, the role of neutrinos as hot dark matter (HDM) particles 



with a large velocity dispersion has been perpetuated in the literature |llMl3| . The approach 
here introduced allows for parameterizing the transition between HDM and cold dark matter 
(CDM) through the analytical expressions for energy density and pressure of a DFG. 

The formalism for a DFG describes a canonical ensemble that, when cooled below a 
critical temperature, presents a condensed phase consisting of quasidegenerate fermion lump 
Ij]- If]]. Self-gravitating fermion systems, for instance, could be conceived through the 
introduction of a distribution function modified with an energy cutoff in phase space, as it 
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happens for a DFG [17H19|. It has also been showed that for a self-gravitating Fermi gas at 
finite temperature, below a critical temperature, it undergoes a first order phase transition 
leading to a compact object with a degenerate core. It is expected that at some stage of 
the evolution of the Universe, primordial density fluctuations have become gravitationally 
unstable forming dense clumps of dark matter that have survived in the form of galactic 
halos [20] . 



The possibilities of forming large structures have stimulated us to investigate the behavior 
of a DFG test fluid in the cosmological background of radiation, matter and ACDM. Due 
to phenomenological reasonable arguments, one should expect a small contribution from 
such test fluid to the cosmic dynamics. In addition, the procedure described in this work 
is complementary to previous studies of massive neutrinos in cosmology, where the study 



of NR and UR fluids (i. e. massive and massless neutrinos) 2lN23| had to be performed 
separately from the beginning of analysis. 

To study the cosmological evolution of a DFG as a test fluid in some dominant cosmolog- 
ical background, our manuscript was organized as follows. In section II we report about the 
usual analytical approach for the Einstein equations for perturbations. Part of our results 
are depicted from section II. We show how to treat analytically the perturbations for rela- 
tivistic and non-relativistic fluids based only on three realistic approximations. In section 
III we adapt such results for obtaining the evolution equations for the density and fluid 
velocity perturbations of a DFG. In section IV we obtain the preliminary results for density 
fluctuation, 5, and fluid velocity divergence, 6, for a DFG test fluid in the RD background 
universe. Extensions to MD and A-dominated (AD) background scenarios are considered. 
It consistently follows analytical continuity conditions. We draw our conclusions in section 
V. 

II. ANALYTICAL APPROACH FOR PERTURBATIONS 

For a homogeneous Friedmann-Robertson- Walker (FRW) flat universe with energy den- 
sity, p(r), and pressure, P(t), one has the following evolution equations for the scale factor, 
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a(r), 



-(i)._|a^ +tf) (2) 

where the dots denote derivatives with respect to the conformal time, r. It follows that the 
expansion factor scales as a oc r in the RD era and as a oc r 2 in the MD era. 

We find that it is more convenient to solve the the first-order perturbed Einstein equations 
in the conformal Newtonian gauge, so that one has 
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where the variables 6* and a are defined as 

(p + p )e = > (p + p )(7 = _ I^) E * , (7) 

and E* = T* — 5jT*/3 denotes the traceless component of T*. 

For a perfect fluid of energy density p and pressure P at a comoving frame, the energy- 
momentum tensor has the form 

TZ = Pg» + (p + P)U»U u , (8) 

where = dx^ /\/—ds 2 is the four- velocity of the fluid. For a fluid moving with a small 
coordinate velocity v i = dx { /dr, v i can be treated as a perturbation like 5p = p — p, 
5P — P — P, <f> and ip. Then to linear order in the perturbations the energy-momentum 
tensor is given by 

T ° = -{p + 5p), 

r = (p + p) Vt = -t; , 

T. = (P + 6P)6) + E* , e; = , (9) 
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where we have allowed an anisotropic shear perturbation S* in T*f. At this point, it is 
pertinent to observe that, since the photons are tightly coupled to the baryons before re- 
combination, the dominant contribution to this shear stress comes from the cosmological 
neutrinos. One notices that, for a fluid, 6 defined in equation (J7J) is simply the divergence 
of the fluid velocity: 6 = i k 3 Vy 

Through the Bianchi identities, the Einstein equation implies the conservation of the 
total energy-momentum tensor. Actually, the energy-momentum tensor of each uncoupled 
fluid is conserved, and obeys the continuity equation in fc-space, 



5 = -(1 + w) (0-30) , (10) 

and the Euler equation, 

n in in 

9 = — (1 - 3w)6 — 9 + — — k 2 5 - k 2 a + k 2 iP, (11) 

a v ' 1 + w 1 + w v ; 

where we have assumed that, for isentropic primordial perturbations, 5P/5p = c\ = P/p = 
w, where c s is the adiabatic sound speed in the fluid. In the non-relativistic fluid description 
for CDM and the baryon components, one sets w ~ 0. The photon and the neutrino compo- 
nents should be more accurately described by their full distribution functions in phase space. 
Herewith we shall develop the procedure for treating perturbations valid for relativistic and 
non-relativistic test fluids in a DFG regime. 

To treat a DFG as a test fluid in some generical cosmological background, we have to 
follow a sequence of three approximations. Let us consider that the phase space distribution 
of the particles gives the number of particles in a differential volume dx 1 dx 2 dx 3 dp 1 dp 2 dp 3 in 
phase space, 

f(p(E, nj ),x\T) dx 1 dx 2 dx 3 dp 1 dp 2 dp 3 /(2ir) 3 = dN. (12) 

where / is a Lorentz scalar that is invariant under canonical transformations. Assuming 
that linear perturbations on the temperature is parameterized by 

T(l + A0rSn 3 -,r)), (13) 

with T = T(t) = T(a), the zeroth-order element of the phase space distribution, 
f{p(E,n j ),x i ,T), is the Fermi-Dirac distribution given by 
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where the factor g is the number of spin degrees of freedom, the Planck and the Boltzmann 
constants, h and k, were set equal to unity, and T = ix {T / /i) denotes the temperature of the 
particles today. In addition, we have introduced into the Eq. (j!4p our first approximation. 
The thermodynamics of our test fluid is constrained by the cosmological behavior for the 
chemical potential, /i ~ /i(a). In the literature, the contribution due to the chemical potential 
fi to the distribution function is often discarded by setting /i = 0. Since we are interested in 
depicting the cosmological behavior of a DFG, we cannot perform the same approximation, 
i. e. we cannot set fi = 0. The simplest way to manipulate it in the Boltzmann equation 
by keeping valid the condition 



dr 



A=0 



da 







(15) 



A=0 



is by assuming that fi/T = fi /T does not depend on a, i. e. that the chemical potential 
follows the same temperature fluctuations on time. Consequently it does not contribute to 
first-order corrections of the Boltzmann equation and the following constraint between / 



and T is maintained 
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?L = E?L = -I ( T ^-\ (16) 
dE p dp E \ dT J ' { ' 

The Fermi momentum of a DFG will thus be written as pp = qp/a and its comoving 
representation as qp = \fp^ — m 2 a 2 , i. e. the ratio g^/m diminishes as the Universe expands. 
If the chemical potential at present, fi , is larger than the particle mass, m, of the test fluid, 
it is possible to depict analytically the transitory regime, from relativistic to non-relativistic, 
for a DFG. 

To go further we must now expand / about its zeroth-order value as 
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(17) 



where the perturbed piece is concerned with our second approximation since it can be written 

as 

T 9f A 1 E 2 df 

f(E,r)dT f(E,r) p dp 

The general expression for the energy-momentum tensor written in terms of the distri- 
bution function and the momentum components is given by 



T£ = J dp4v.dp z {-Det\g^\)- xl2V ^fW,n 3 \x\T\ 



(19) 



In the conformal Newtonian gauge, (— -Detfg^j,]) -1 / 2 = a~ 4 (l — ip + 30), dp 1 dp 2 dp 3 = (1 — 
3(p)q 2 dq dfl, and p = — (1 + ip)e. Keeping e as a(r) times the proper energy measured by a 
comoving observer and introducing the tilde-vaxiables as 



(27T) 



(20) 



to linear order in the perturbations one writes 

T ° =- / ^ 2 £/(£,t)(1 + M/), 



rj-ii 



dpp 3 f(E,r)ty n j 



dp^f{E, r ){l + **)n%. 



(21) 



Noticing that the angular dependence in terms of x i and n can be decoupled from the q- 
integration through Eq. (I18p. the tt/de-perturbations to the energy-momentum tensor can 
be written as 
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A(x\n j; T), 

A(x i ,n j ,T) Tlj, 

A(x\n j ,T)n t n :j 



(22) 



One can notice that we have eliminated the explicit dependence on the metric pertur- 
bations. Rewriting the integrands in terms of the (pseudo) comoving energy, e = a E = 
a (p 2 + m 2 ) 1 / 2 = (q 2 + c^m 2 ) 1 / 2 , and of the comoving momentum, q = pa, and integrating 
by parts, it results in the following expressions for the t«/de-perturbations to the energy- 
momentum tensor 
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5T 



a(r)" 4 
a(r)- 4 
a(r)" 4 



m 2 a(r) 2 



ef(E 7 r) + A dqq 2 ef(E,r) 



dqq 3 f(E,r) + 2m 2 a(r) 



A(x i ,n j ,T), 
A(x i ,n j ,T) Tlj, 



dqq 2 ef(E,T)+ / dq — f(E, r) 



A(x i ,n j ,T) njij 



where it is easy to identify the average energy density and pressure respectively as 

1 



P(r) 
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dqq 2 ef(E, T ), 



(23) 



(24) 



and 
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P(t) = / dq — f(E,T), (25) 

V ; 2vr 2 a(r) 4 J q 3e JK h K J 

and where A(a; i ,%,r) reads J\f(k i ,n j ,T) in the Fourier space, i. e. 

A(s\nj,T) — t — r / dA; 3 exp [ik,x'] J\f{k\h h T). (26) 
(2tt) 3 J 

When particles are ultra-relativistic, i. e. e = q, the above results reduce to that obtained 



for massless particles |21|, |22j. Differently from that observed when one sets the ultra- 
relativistic approximation on f(E,r), the spectrum does not remain Planckian. However, 
under certain particular circumstances, the local dependence on e(q) and q remains easy to 
treat analytically. The g-integration can be performed separately and the dependence on x\ 
n and r is kept active. 

Now, the dimensionality of the problem can be reduced by noticing that the evolution 
depends only on the direction n through the angle related to k • n — cos <p, which is a natural 
consequence of the isotropy of the homogeneous background. We perform the P s -Legendre 
expansion with respect to cos(y?), 

oo 

Af(x\h J7 T) = ^(-l) s (2s + l)iV s (fc\T)P s (cos^), (27) 

that results in 

N = -L J dQ Af, N 1 = ^- J dfl cos ipAf and iV 2 = - A J dQ ^cos 2 <p - ^ J Af, (28) 

for the first three n-poles. 

Introducing the above results after performing the angular integrations over the tilde 
variables, the Eqs. (123|) result in 



<5T ° Sp = N (a P (t) + £^- 2 J dq e f(E, r^j 



dqq 3 f{E,r) + 2m 2 a 2 dqqf{E,r) 



k N 

ik m 5V m =► (p( T )+P(r))e = ^ 
- (l/3)^)S* (p(r) + P(r))a = 2N 2 (p(r) + P(r)). (29) 

The results of this section allows one to perform preliminary analytical computations 
for the evolution of perturbations for several classes of relativistic and non-relativistic test 
fluids. As an example, in the section IV we shall obtain such results for a DFG in the RD 
universe and some extensions to MD and AD scenarios. 
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III. PERTURBATIONS AND EINSTEIN EQUATIONS FOR A DFG 



The third approximation is the simplest one. It sets the DFG cut-off on the distribution 
function f(E, r), which results in obtaining analytical expressions for the elements described 



in Eq. 



In the literature 



21 



22 



24|, however, one can notice that f(E, r) is often 



approximated by f(p) = (exp [p/T] + Besides not resulting in the expected smooth 

transition between relativistic to non-relativistic regimes, it is out of the scope of the study 
of a DFG cosmological test fluid. From such a point, the perturbations to a DFG can be 
analytically expressed without additional assumptions over f(E, r). 

In the limit where T tends to in Eq. (fTTj) . the Fermi distribution f(E, r) becomes 
a step function that yields an elementary integral with the upper limit equal to the (co- 
moving) Fermi momentum, qp- Introducing the variable y(a) = A/a, with A = q F /m = 
(q F /T )/(m/T ) to parameterize the kinetic energy with respect to the rest mass, the equa- 
tion of state of the DFG can be expressed in terms of elementary functions of y = y(a(r)) 



as 



p(y) 
P{v) 

n(y) 



g m 

16 7T 2 

g m A 

16 7T 2 



y(V + 1) Vv 2 + 1 - arc sinh (y) 



y{^y 2 - l)vV + 1 + arc sinh (y) 



q m , 

6 7T 2 



(30) 



from which ultra-relativistic and non-relativistic regimes can thus be analytically connected 
during the cosmological evolution. 

The pressure and energy density are then related by 

dp(y) p(y) + P(y) 



vV + i- 



(31) 



d(mn(y)) rnn{y) 

From Eqs. ( 1301) . the relations of the density fluctuation, 5, the fluid velocity divergence, 
8, and the shear stress, a, with the multipoles of lower order, iV , N± and iV 2 can be then be 
rewritten as 



m 



5p = N y(l + yT 2 — , 

7T 2 



9 = 3k N, 



VT+Y 2 
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2N 9 



(32) 
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as well as 



5P = N y 3 ^T+Y 2 



m 
3tt~ 2 



(33) 



In Fig{T]we describe the most relevant perturbation variables for a DFG, by quantifying 
their relative magnitudes with respect to the energy density fluctuation, and by comparing 
it with the results for radiation perturbation. 

By introducing the above simplifications given by Eqs. fj30|) into Eqs. (llOflip for the 
evolution of perturbations one thus obtains 
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dr 
and 
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dr V P(y(r)) 



P{y{r)) 



p(y(r)) + P(y(r)) 



p(y(r)) + P(y(r)) dr dr^ 



k 2 5 + k 2 0, (34) 



(35) 



through which, under simplified assumptions for the background cosmological scenarios that 
set the dynamical behavior of the Hubble parameter, H(a), we obtain the DFG test fluid 
evolution. 

One should obtain the same results departing from the Boltzmann equation in fc-space, 



— + i-(k-n)y + — 

or e dmq 
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c 



once the condition given by ffl5l) had been set. For instance, the equation for the shear 
stress, a = 2N 2 stays as 

iL(£ iVl ) +fc £ [(T _ J V„] = 0. (37) 



IV. DFG TEST FLUID PERTURBATIONS 



Deep inside the RD era, one sees that the radiation component can be approximate d by 



a perfect fluid of sound speed c s , due to the tight coupling between baryon and photons [24]. 
This self-gravitating fluid oscillates inside the sound horizon. The cold dark matter evolves 
as a test fluid in this background, until the time (close to equality) at which its gravitational 
back-reaction becomes important. In this limit, and in the absence of anisotropic stress 
(0 = ip), it is possible to derive an analytic expression for the evolution of the metric 
perturbation <p 1 . The full continuity, Euler and Einstein equations for the RD background 
results in the following expression for the scalar perturbation (see, for instance, [21 . H]) 



[OJT 
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+ sin wr| 



(38) 
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where lo is proportional to the Fourier wavenumber k, u — k/y3, and r is the conformal 



time related with the proper time through the scale factor by dr = dt/a(t) 22J. We just 
have noticed that the simplest way to describe the subsequent results is in terms of the 
dimensionless parameter ut, i. e. the conformal time weighed by u. 

The solution proportional to C represents the growing adiabatic mode, the other one pro- 
portional to D is the decaying adiabatic mode, which is suppressed during the cosmological 
evolution. 

By substituting the above expression for 7 into the Eqs. ( )34|) - (l35|) one can obtain 5\ and 
9\ numerically for a DFG test fluid. For preliminary results, we have approximated the Uni- 
verse composition by VL 1 = 0.01, fl m = 0.29, and — 0.70. In order to perform a generical 
analysis in terms of non-relativistic to ultra-relativistic regimes (A) and cosmological scales 
(ur), we have considered different values for the parameter A in resolving Eqs. (j34p - (!3"5|) . 
The results can be depicted from Figs. |2] and [3] where the dependence on the relativistic 
regime through A does not obey a regular behavior, i. e. the most subtle effect on the DFG 
cosmic fluid is the suppression of the growing modes while the RD universe is cooling. In 
some sense it is unnatural since one should expect DFG approximating the non-relativistic 
matter behavior (growing modes) as the universe cools. We compare the results to those for 
matter (blue line) and radiation (red line) obtained from direct resolution of Eqs. ( jTOlfTTT) . 
respectively in the NR and UR limits, 

L. = - ( e. 



and 



e m = --9 m + 3co 2 (P, (39) 
a 



= (fcy + 40 



3 

6> 7 = ^ 2 <5 7 + 3w 2 0. (40) 

Turning back to the relation between the conformal time, r, and the scale factor, a(r), for 
the RD era, one has 

a 7 (r)~ft 7 / 2 r, (41) 
and, by adopting adequate continuity conditions, one sets 

(42) 
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and 

1 ( 3 n 1 ' 2 y 1 

a A (r) r — / r ~ tt- - r ) , 43 

respectively for the MD and AD eras. One thus finds the conformal time for the equality 
between radiation and matter densities, r eq , for the equality between matter and cosmo logical 
constant A densities, r m A, and for the present A dominance at a = 1, tao, 



t, 



1/2 



2 

TmA ~ i/3 i/ 6 

3 1 1 

TAO ~ ^T m + -T eq - —j- 2 . (44) 

In addition, the analytical solutions for the metric perturbations <p(r) related to the MD 
and AD regimes, are respectively given by 

m (wr) = C m {uj) + m<y J , (45) 
(ur)° 

and 

<j) m {ujT) = C A M(a A (r))- 3 + D A (uj)(a A (r))-\ (46) 

By fine-tuning the transition between RD, MD and AD solutions through adequate conti- 
nuity conditions between the growing adiabatic mode of 7 and the complete solution for 
<f) m and between their first derivatives, both at T eq , and subsequently doing the same for the 
complete solutions of m and 0a at r m , one can easily find the coefficients C m , Ca, D m , and 
D A . 

This corresponds to the simplest strategy for obtaining the first approximation for a 
transference function between radiation, matter and A eras, since it allows one to build 
continuous solutions for the density fluctuation 5\ and fluid velocity divergence 9\ of the 
DFG test fluid. In the Fig. H] we have illustrated the results for 5\, in correspondence with 
the A values introduced in the previous results for classifying the DFG test fluid according 
to its relativistic regime. The cut-off in each plot is in agreement with the values considered 
in (ED. 

One has to notice that the relation between 9 and u are relevant for determining the be- 
havior of the density fluctuations. The first stage (i. e. at RD era) of growing perturbations 
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depends on how smaller is A, i. e. they are favored by non-relativistic configuration regimes 
of the DFG test fluid. Otherwise it also depends on scales parameterized by the weight 
factor u in the x-axis for the conformal time. When the perturbation scale is well inside 
the Hubble horizon, the oscillating behavior (RD) followed by a strong suppression at the 
last stage (AD) is favored. For large scales, the growing perturbations may or may not be 
suppressed at the last stage. Such a qualitative preliminary analysis allows for reheating the 
discussion on the possibility of forming DFG structures at different cosmological eras, and 
deserve some attention, since the behavior of DFG perturbations does not mimics matter 
perturbations exactly, with the exception of the totally NR limit. 

V. CONCLUSIONS 

An analytical procedure for obtaining the evolution of perturbations for a class of rela- 
tivists and non-relativistic test fluids was discussed. The density fluctuation, 5, the fluid 
velocity divergence, 9, and an explicit expression for the dynamics of the shear stress, a, 
were obtained for a DFG in the background regime of radiation, and some complementary 
results were obtained in the case matter domination and ACDM cosmological eras. The 
comparisons with previously known results from radiation and matter perturbations are im- 
mediate from our analysis and show the consistence with the DFG approach in the sense 
that an equivalent equation of state can smoothly transit from the UR to the NR limits. In 
particular, one can observe that a DFG, even in the non-relativistic regime at early times 
preconize the growing of perturbations in the same intensity as matter. However, for inter- 
mediate regimes, the growing adiabatic modes are relatively suppressed after the beginning, 
turning to an oscillatory behaviour. Also for the fluid velocity divergence, the transition 
from NR to UR regimes can be continuously observed. In this sense, our results become a 
convenient tool for performing preliminary tests for HDM to CDM test fluid configurations. 
Moreover, the DFG approach can be used to set a correspondence between the averaged 
neutrino mass at present (a = 1), thorough the NR limit of the energy density, and the 
Fermi momentum, pp, through the UR neutrino energy density at the RD era. 

Concerning the neutrino degeneracy, one can ascertain that the value of the chemical 
potential of electron neutrinos (/i e /T ) should be close to —1 in order to explain with an 
acceptable confidence level the phenomenology related to the observed variation of deuterium 
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by roughly an order of magnitude [23]. It would induce a variation in total energy density 
during the RD stage, which is excluded by the smoothness of CMB. However, this objection 
could be avoided if there was a coincidence between different leptonic chemical potentials 
such that in different spatial regions they had the same values but with interchange of 
electronic, muonic and/or tauonic chemical potentials [30]. For further details and discussion 



of energy dependence and effective chemical potential one might address the papers [30l. |31| 
where the main point is that they are not obtrusive to our results. 

It should also be mentioned that cosmological models with a large excess of neutrinos - or 

n 

antineutrinos - for example, have been discussed [25] in this context. In such models there 



is a large neutrino background behaving like a DFG. The equilibrium distributions show 
that for fi u 7^ 0, and \/jl v \ 3> kT, there is an excess of particles. During the RD phase of the 
early universe, a non-zero chemical potential for neutrinos, fi u , can give rise to an excess of 
particles over antiparticles or vice-versa. The massless neutrinos (or antineutrinos) would 
then form a Fermi sea with degenerate particles filling the Fermi levels up to pf- However, 
the exact neutrino number density is not known as it depends on the unknown chemical 
potentials for the three flavors. (In addition there could be a population of sterile neutrinos 
29I ]. hypothesis that we did not discuss here.) In this context it is also frequent to notice the 
term neutrino stars as a generic name for any degenerate fermion star composed of neutral 

ml 
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weakly interacting fermions, e.g., neutrinos or supersymmetric fermionic partners 

Finally, the hypothesis of a small fraction of dark matter components evolving cosmologi- 
cally as a DFG test fluid in some dominant cosmological background can also be speculated. 
The idea is based on the fact that for any type of fermionic dark matter, the average phase- 
space density cannot exceed the phase-space density of the DFG. In this case, instead of 
specifically considering structures like exotic stars or exotic planets, one should consider the 
possibility of forming extended overdense regions. This study can be relevant for massive 
neutrinos in dark matter models and for collisionless self-gravitating systems experiencing 
some kind of relaxation. Our results could thus be extended to support the discussion of the 
formation of the halo of our galaxy assuming that it consists of a gas of massive degenerate 
fermions in hydrostatic and thermal equilibrium at finite temperatures. Such a formulation 
is consistent with the precursive idea of a self-gravitating Fermi gas model introduced to 



explain the puzzling nature of white dwarf stars 



32 



331 ] . This fermion gas model is expected 



to imply the existence of supermassive compact dark objects at the galactic center. 
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FIG. 1: Evolution of the DFG relative isotropic pressure perturbation dP/p, fluid velocity diver- 
gence, 9, and anisotropic stress, a, in units of the density perturbation, 5, in dependence on the 
inverse of the scale factor, a -1 , weighed by the DFG factor A. In the non-relativistic limit, i. e. 
when A <C 1, one recovers the values for matter perturbation. In the ultra-relativistic limit, i. 
e. when A> 1, one recovers the values for radiation perturbation. The comparison to radiation 
perturbations follows in the second plot. 
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FIG. 2: Evolution of radiation (<5 7 ), matter (5 m ) and DFG (5\) energy density fluctuations in 
dependence on the weighed conformal time, cot, during the RD regime. The parameter A summa- 
rizes the dependence on the relativistic regime and on the cosmological scale. By varying A one 
describes how relativistic is the regime; and by varying to one sets the relation to the cosmological 
scale. The results are for A = 0.01, 0.1, 1 and 10. The initial condition was set arbitrary since we 
are depicting only relative values. 
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FIG. 3: The fluid velocity divergence, 9, in dependence on the weighed conformal time, wr, during 
the RD regime. The results are also for A = 0.01, 0.1, 1 and 10. 
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FIG. 4: The density energy fluctuations, 5\, in dependence on the weighed conformal time, ur, 
respectively for RD, MD and AD regimes. Time is evolving from the beginning of the RD era. The 
values for the parameter A are in correspondence with the previous figures. 
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